#include <stdio.h>
#include <stdlib.h> 
#include <math.h>
#include "system.h"

// read in system information
void Readdat(void)
{ 
    NumberOfParticles=300;
    NumberOfSteps=6000;
    Temperature=1.428;
    Tstep=0.001;
    Density=0.728;
    NumberOfInitializationSteps=2000;
    Box=pow(NumberOfParticles/Density,1.0/3.0); /* Box Length */
    
  if(NumberOfParticles>Maxpart)
  {
    printf("Maximum number of particles is : %d\n",Maxpart);
    exit(1);
  }
 
  // Calculate Some Parameters
  CutOff=2.5;
  Ecut=4.0*(pow(SQR(CutOff),-6.0)-pow(SQR(CutOff),-3.0));
 
  // print information to the screen
  printf("Molecular Dynamics Program\n");
  printf("\n");
  printf("Number of particles   : %d\n",NumberOfParticles);
  printf("Boxlength             : %f\n",Box);
  printf("Density               : %f\n",Density);
  printf("Temperature           : %f\n",Temperature);
  printf("Cut-Off radius        : %f\n",CutOff);
  printf("Cut-Off energy        : %f\n",Ecut);
  printf("Number of steps       : %d\n",NumberOfSteps);
  printf("Number of init steps  : %d\n",NumberOfInitializationSteps);
  printf("Timestep              : %f\n",Tstep);
 
}
